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We study the metal-insulator transition in a tight-binding one-dimensional (ID) model with long- 
range correlated disorder. In the case of diagonal disorder with site energy within [— 2.] and 
having a power-law spectral density S(k) oc k~ a , we investigate the competition between the disorder 
and correlation. Using the transfer-matrix method and finite-size scaling analysis, we find out that 
there is a finite range of extended eigenstates for a > 2, and the mobility edges are at ±E C = 
±|2 — W/2\. Furthermore, we find the critical exponent v of localization length (£ ~ \E — E c \~ v ) to 
be i/ = l + lAe 2 ~ a . Thus our results indicate that the disorder strength W determines the mobility 
edges and the degree of correlation a determines the critical exponents. 

PACS numbers: 72.15.Rn, 71.30.+h, 73.20.Jc 



I. INTRODUCTION 

Anderson's localization theory points out that in a sys- 
tem with uncorrelated diagonal disorder, all one-electron 
states are spatially localized with an exponentially de- 
caying envelope with a characteristic localization length 
£, when the disorder strength is larger than a critical 
valued From the scaling theory^ it is also well-known 
that an infinitesimal disorder can cause localization of all 
states in one and two dimensional systems in the ther- 
modynamic limit. In recent years, it was found that spa- 
tial correlation of disorder potential played an important 
role in the nature of long-range charge transport of low- 
dimensional systems. In the presence of short-range cor- 
related disorder, there exists a set of discrete resonant 
energy levels of extended states i&^&^ i 8 ' 9 ' 10 ! 11 ' 12 For ex- 
ample, a single extended eigenstate was found in the ran- 
dom dimer mo del j_L2_12i which was verified by measuring 
microwave transmission on the semiconductor supcrlat- 
tices with intentional correlated disorder* 1 ^ These models 
with short-range correlated disorder do not possess a true 
disorder induced metal-insulator transition in the ther- 
modynamic limit. More recently, a one-dimensional (ID) 
disordered model with long-range correlation considered 
by de Moura and Lyra has arisen a great interest, as it 
can result in a continue band of extended states under ap- 
propriate conditions. 14 The transition from localization 
to derealization occurred at the critical points was ex- 
amined experimentally in a single-mode waveguide with 
inserted correlated scatters* 1 ^ More recent experiments 
in ultra-cold atoms/Bose-Einstein condensate push the 
research in this direction further^ 

In the disordered systems, the critical parameters, such 
as the critical disorder strength W c , mobility edge E c , 
critical exponents i>e, vw (defined by the dependence of 
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localization length on W or E near the critical points 
£ ~ \W - W e \ vw or £ - \E - E C \ VE ) are of partic- 
ular importance and interest, since they determine the 
nature of the localization-delocalization transitions and 
the universal properties of disordered systems. Even for 
three-dimensional (3D) Anderson model which has been 
studied for more than forty years, the accurate deter- 
mination of the critical parameters is still under recent 
investigations^ where the critical parameters for differ- 
ent lattices are calculated and it was found that ve is 
close to vy/. While for 3D systems with scale- free disor- 
der, it was found that ve ^ vw and the extended Harris 
criterion is obeyed. — 

In spite of many studie^^^^^ 2 .^^^^^^^ 
on the ID systems with correlated disorder, many of 
them are in a perturbative level and there are still 
several important questions to be answered, such as 
the positions of the mobility edges and the value of 
the critical exponent v. The determination of these 
critical parameters and the comparison with the results 
of 3D models mentioned above will definitely deepen our 
understanding on disordered systems, especially the de- 
pendence on dimensionality and correlation. Moreover, 
the correlated disordered system is a kind of system 
which is intermediate between ordered system and pure 
disordered system. The investigation of the universal 
properties of this type system and comparison with 
quasi-periodic system is of fundamental importance. In 
addition, the accurate values of these critical parameters 
(mobility edges) have important applications in charge 
transport in low dimensional (correlated) disordered 
systems, such DNA molecules* 2 ^, Bose-Einstein conden- 
sate in 1-D optical lattices with speckle potential^. In 
the present paper, we calculate these critical parameters 
for the ID systems with the correlated diagonal disor- 
der. And we also resolve some inconsistencies in the 
literatures* 1 ^ One of the advantages of our study is that 
we are able to investigate localization-delocalization 
transition and obtain the accurate critical parameters 
in ID system with large size and the boundary effects 



2 



is negligible, unlike the 3D system where all the sizes 
in three dimensions can not be very large due to the 
computational limitation. 

We first investigate the localization properties of the 
ID tight-binding model with long-range correlated diag- 
onal disorder (the hopping constant is set to be unity). 
The diagonal on-site energies {e^} are distributed in 
[-W/2, W/2] with W being the strength of disorder. The 
site energies have a power-law spectral density S(k) oc 
k~ a with a > 0. The function S(k) is the Fourier trans- 
form of the two-point correlation function (cjCj) and k 
is related to the wavelength A of the undulations on the 
random parameter landscape by k — We find that 

the critical disorder width W c = 4 — 2\E\ at fixed energy 
E for a > 2. Moreover, we can determine the positions 
of the mobility edges ±E C = ±|2 — ^-|, which separate 
the extended and localized energy eigenstates. The ef- 
fective energy band width B of the extended states has 
the linear relationship B = 4 — W ioi a > 2. Thus, a 
phase diagram in the (E, W) space is obtained. We also 
discuss the a-dependence of the critical exponent v of 
the localization length of the eigenstates. Interestingly, 
we find that the disorder strength W determines the mo- 
bility edges and the degree of correlation a determines 
the critical exponents. Our results also indicate that the 
nature of a correlated disordered system is somehow be- 
tween that of the pure random (without any correlation) 
and pseudo-random (quasi-periodic) systems. Our ap- 
proach is non-perturbative and our results are hard to 
be obtained by perturbative calculations. 

The organization of the paper is as follows. In the 
next section, we introduce the ID disordered systems 
with long-range correlation and the basic approach in 
our calculation. In section III, we present our results of 
the mobility edge and critical exponents. The paper is 
summarized in section IV. 



II. MODEL AND THE BASIC APPROACH 

We consider a ID tight-binding Hamiltonian with long- 
range correlated disorder 

H = J2^\n)(n\ + J2 t i\n}(n + 1| + |»+l)<n|), (1) 

n n 

where \n) is the Wannier state localized at site n with 
on-site energy e„, and t is the nearest- neighbor hopping 
amplitude. We set t = 1 for simplicity. The on-site 
energies {e n } can be constructed by the Fourier filter- 
ing metho d 24 i 25 ' 26 as follows: (i) generate a random se- 
quence {u n } with a Gaussian distribution; (ii) get its 
Fourier components {uk} using the fast-Fourier trans- 
formation method; (iii) establish a new sequence {e^.} 
by the relation tk = k~ a l 2 Uk] (iv) construct the se- 
quence {e n }, which is the inverse Fourier transform of 
{efc}; (v) adjust the scale of the sequence {e„} reach- 
ing to [— W/2, W/2}. It can be checked that the disor- 
dered on-site energy {e n } is long-range correlated with 



the power spectrum S(k) oc J2 m (e n £n+m)e' lmk oc k~ a . 27 
The exponent a characterizing the correlation reflects the 
roughness of the energy sequence. The larger the value 
of a is, the smoother the energy landscape is. 

The Schrodinger equation for the wavefunction ampli- 
tudes 1p n is 

e n 4>n + 1pn-l + 1pn+l = Elp n . (2) 

Here, E is the eigenenergy. Our interest focuses on the 
critical behavior in the thermodynamical limit. A useful 
way to find the critical points is based on the finite size 
analysis of the normalized localization length A = £/N 
If A increases as the system size N — ► oo, which im- 
plies that £ would exceed N when N is large enough, 
the electron stays on a extended state; Otherwise, the 
state is localized, because A is a decreasing function of 
N and smaller than unity in the thermodynamic limit. 
So the N dependent of A = £/N can give the information 
of the localization-delocalization transition. Thereby, we 
can use this quantity A to determine the critical disorder 
strength at fixed energy and the critical energies at fixed 
disorder strength. 

In the regime of weak disorder, it was founcU^ that the 
localization length £ for a state of energy E = 2 cos /i 
(not near band center and band edge) is determined by 
the correlation function of the disorder potential 

r 1 = two*), (3) 

where 70 = 8sin a (/J ) , = 1 + 2 Em=i C0S ( 2 H9m, 

(e n e™+m) = £o<7m, (4) = e o- For our disordered sys- 
tem with long range correlation with the spectral density 
S(k) oc k~ a , the localization length has the normalized 
form A = £/N - N a ~ 2 , for 1 < a < 2; A ~ 1/JV, for 
a < 1 . Thus we see that in the weak disorder regime the 
states are localized for a < 2 and a = 2 is a critical value 
for the appearance of extended states. 

To determine the full phase diagram and study the crit- 
ical parameters, we use the transfer matrix method. Us- 
ing the two-component vector <!>„ = (tp n , i/) n - i) T , where 
the superscript T denotes the transpose, Eq. (2) can be 
written in a recursive form 

(V )-*•(&) (4) 

with 

«>=(V' "o 1 )' (5) 

We obtain the transfer-matrix equation &n+i = TW^i, 
with Tjv = YiiLi The localization length £ at energy 
E is the inverse of the Lyapunov coefficient (7 > 0), 
which is the largest eigenvalue of the limiting matrix 
limAr^oo ln(T N Tj f ) 1 / 2N 2& The large numbers appeared 
in the calculations have been took care of by dividing 
a large number and multiplying it at the end of calcu- 
lation. The reorthogonalization method 29 has also been 
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FIG. 1: The normalized localization length A = £/N as a 
function of the system size iV for fixed E = 0. a is taken as 
3.5 and 1.5 in (a) and (b), respectively. 



used and the same critical parameters are obtained. All 
the values of 7 in this paper are based on a geometrical 
mean of 10 4 disorder configurations. 



III. RESULTS AND DISCUSSIONS 
A. Critical disorder strength and mobility edges 

We perform the calculations of the normalized localiza- 
tion length A for different values of a and energies. Fig- 
ure 1 shows plots of the normalized localization length 
A(W) versus N for E = 0, a = 1.5 and a = 3.5, re- 
spectively. The system size N ranges from 2 13 to 2 17 . 
We notice that A is monotonically decreasing with N for 
a = 1.5; while the picture is different for a — 3.5: the 
values of A rise with increasing N for W < 4, whereas 
the values of A decline with increasing N for W > 4. 
The size independence of A at W — 4 indicates a critical 
point of a continuous phase transition. 

In Fig. 2, we display the dependence of A on the disor- 
der strength W for the state of E = 1 with typical values 
of a. One may observe that there emerges an intersection 
point for different system size when a — 2.5, which cor- 
responds to the critical point. This critical point doesn't 
exist for a = 1.5. From the inset of Fig. 2, we could ex- 
plicitly note the different dependence of A on the system 
size N in the vicinity of W c = 2. It tells us the state 
with E = 1 is localized for a = 1.5 and is extended for 
a = 2.5 and W < 2. 

Table 1(a) lists the critical values of W c for various 
energies. It can be seen that the critical disorder strength 
is independent of a whenever a > 2. At the same time, 
we find a simple relationship between W c and E as: 



2\E\. 



(6) 



For the special state at band center with energy (E = 0), 
the critical value W c = 4, equals to the bandwidth. This 
conclusion for E = state was also reached by H. Shima 
et alZL 

Using the same method, we also obtain the critical en- 
ergies ±E C at fixed W, as shown in Fig. 3. The results are 
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FIG. 2: The normalized localization length A as a function of 
the disorder width W at fixed E = 1 for systems with different 
size, a = 1.5 and a — 2.5. Insert: The detail behavior of 
A(VF) near the critical point W c ~ 2. 




FIG. 3: The normalized localization length A as a function of 
E at fixed W = 2 for a — 3.0. Detail behavior of A-(E) near 
the critical point is shown in the inset. 



given in Table 1(b) for a > 2. The critical energies/the 
mobility edges can be determined by the following equa- 
tions: 



±E r = ±12 



w 



I («>2). 



(7) 



The eigenstates between — E c and E c are delocalized. 
Therefore, the effective bandwidth B of the delocalized 
states is (for a > 2) 



B =\E c -(-E c )\=4-W (a>2). 



(8) 



Figure 4 shows the phase diagrams in the (W, E) and 
(E, a) plane. The critical value a — 2 is also consis- 
tent with the perturbative calculation discussed above. 
We emphasize that mobility edges or the effective band- 
width is only dependent on the disorder strength W as 
long as a > 2. Also, in the limit of W — 0, the ef- 
fective bandwidth B becomes 4 which is the well-known 
result in the ID ordered crystal lattice. It is interesting 
to compare with the pseudo random systems with quasi- 
periodic on-site energy e n = \ cos(nQn a ) € [— -y, ^f], 
Q an irrational number, < a < 1. In this model, the 
critical value W c = 4 — 2\E\ and the mobility is also 
at ±|2 — W/2\, independent of The linear relations 
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FIG. 4: Phases diagrams in the (W, E) and (E, a) plane. The 
delocalized and localized states are separated by W = 4 — 2\E\ 
and a = 2. Here t = 1. We take = 1 in (E, a) plane as an 
example. 
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FIG. 5: Scaling function of In A for a = 2.5. Data from 
various chains, whose size N ranges from 2 13 to 2 17 , overlap 
each other. As an example, we give the optimized values of 
E c and v in the case of W = 2.0. 

between the mobility edge E c (critical disorder strength 
W c ) and disorder strength W (energy E) are simple and 
neat, further more they have a direct consequence on the 
critical exponents as will be discussed in the next sub- 
section. 



B. Critical exponents based on finite-size Scaling 
Analysis 

In the light of the finite-size scaling analysis, the nor- 
malized localization length A for finite system may be 
related to the localization length ^ for infinite system 
by the following scaling lawj2& 

InA = /A, (9) 

soc 

where the parameter £ x varies as ^ = \W — W c \~ Vw 
or Coo = \E T E c \ ~ VB in the vicinity of the critical point. 
Unlike three dimensional systems with usually small val- 



ues of scaling variable (i.e. the size of cross section), 
in our one-dimensional systems, the irrelevant scaling 
exponents^ are of little effect due to scaling variable with 
very large values, i.e., the size of the system. 

Firstly, we expand the scaling function in polynomial 
form: 

In A = cio + a\X + &2£ 2 + • ■ • + a n x n . (10) 

Here x defined as i e \W - W c \~N xlvw or x = \E =F 
E c \N 1 l l ' E , is a nonlinear combination of W(E) and N 
with the parameters W C (±E C ) and v to be determined. 
Secondly, based on Eq. (10), we can obtain the values 
of W C (±E C ) and v by fitting A which we got from previ- 
ous calculation. The expansion in the scaling function 
must be carried up to n, that makes sure the degree 
of confidence reaches to 99%. As shown in Fig. 5, we 
have an excellent scaling curve, accompanied by the op- 
timized values of W C {±E C ) and v. The critical values for 
W C (±E C ) obtained from finite-size scaling analysis are 
identical to those obtained in section III (A). The values 
of critical exponent v and their errors estimated from the 
fitting procedure are given in Table II. It is clear that the 
critical exponent v is independent on W c or ±E C . With 
increasing the correlation (a), v decreases to 1. One may 
note that ve = vw- Actually it is a direct consequence 
of the linear relation between the disorder strength and 
mobility edges. From £ ~ \W - W c \~ vw and Eqs. (6) 
and (7), one finds that £ ~ \E - E c \~ Vw . Thus one has 
ve = v\y . It is interesting to compare our results with 
those in 3D systems. For 3D Anderson model (without 
correlation) vw is close ve, 17 while in the 3D system 
with scale free disorder, vy/ ^ v E - — 

Then wc find the dependence of v on a follows the rule: 

v= l + ^iexp(^ 2 (2-a)), (11) 

with A\ is approximately 1.4 and — 1 as seen in Fig. 
6. We have found that v = 1 for infinite large correlation 
of a. It is interesting to note that the case with infi- 
nite large long-range correlation (a — > oo) looks like the 
quasiperiodic system, since the position of the mobility 
edges and the critical exponents (1) in both systems are 
the same. — In another limit situation a = 2, v = 2.4, 
which is not easy to be obtained from direct calculation 
due to the critical slowing down phenomena. This result 
is apparently different from v = 2 at a = 2 mentioned 
in Ref. — . Our results are based on finite size scaling. 
The data show very good scaling behavior as shown in 
Fig. 5. The critical exponents and mobility edges are ob- 
tained simultaneously by using fitting formula Eq. (10). 
The mobility edges obtained in this way are identical to 
those obtained in section III(A). Therefore our results are 
reliable. It is interesting that the mobility edges are de- 
termined only by the disorder strength W and the critical 
exponents are determined only by the degree of correla- 
tion a. 
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1(a) Critical values of disorder strength W c for different a and energies E 





E=-1.5 


E=-1.0 


E=-0.5 


E=0.0 


E=0.5 


E=1.0 


E=1.5 


a 


w c 


w c 


w c 


w c 


w c 


w c 


w c 


2.50 


1.000(1) 


2.001(1) 


3.001(2) 


4.001(2) 


3.001(2) 


2.001(2) 


1.000(2) 


3.00 


1.0000(4) 


2.0000(5) 


3.0001(5) 


4.0001(6) 


3.0001(6) 


2.0001(6) 


1.0001(4) 


3.50 


1.0000(1) 


2.0000(1) 


3.0000(1) 


4.0000(1) 


3.0000(1) 


2.0000(1) 


1.0000(1) 


4.00 


1.0000(1) 


2.0000(1) 


3.0000(1) 


4.0000(1) 


3.0000(1) 


2.0000(1) 


1.0000(1) 



1(b) Critical values of energies ±E C for different a and disorder strength W 



a 


W=3.0 


W=2.0 


W=1.0 


-E c 


E c 


-E c 


E c 


-E c 


E c 


2.50 


-0.5004(6) 


0.5003(5) 


-1.0001(4) 


1.0002(5) 


-1.5001(5) 


1.5001(5) 


3.00 


-0.5000(2) 


0.5001(2) 


-1.0000(3) 


1.0000(2) 


-1.5000(1) 


1.5001(2) 


3.50 


-0.5000(1) 


0.5000(1) 


-1.0000(1) 


1.0000(1) 


-1.5000(1) 


1.5000(1) 


4.00 


-0.5000(1) 


0.5000(1) 


-1.0000(1) 


1.0000(1) 


-1.5000(1) 


1.5000(1) 



TABLE I: The values of critical disorder strength Wc and critical energies (mobility edges) ±E C obtained by the coordinates 
of the intersection points in the cases a > 2. The errors for the last digits are listed in the parentheses. 



11(a) Critical values for disorder strength W c and exponent vw for different a and energy E 



a 


E=0.0 


E=0.5 


E=1.0 


E=1.5 


W c 


v w 


Wc 


vw 


Wc 


vw 


W c 


vw 


2.25 


4.0015(5) 


2.15(5) 


3.0013(3) 


2.13(4) 


2.0008(7) 


2.12(6) 


1.0002(6) 


2.12(6) 


2.50 


4.0007(5) 


1.86(4) 


3.0007(3) 


1.83(3) 


2.0004(3) 


1.84(3) 


1.0003(3) 


1.81(3) 


2.75 


4.0003(3) 


1.65(2) 


3.0003(2) 


1.65(2) 


2.0001(2) 


1.66(2) 


1.0000(2) 


1.66(3) 


3.00 


4.0002(1) 


1.51(1) 


3.0001(1) 


1.50(1) 


2.0001(1) 


1.50(1) 


1.0000(1) 


1.50(2) 


3.25 


4.0000(1) 


1.40(1) 


3.0000(1) 


1.40(1) 


2.0000(1) 


1.40(1) 


1.0000(1) 


1.41(2) 


3.50 


4.0000(1) 


1.31(1) 


3.0000(1) 


1.30(1) 


2.0000(1) 


1.31(1) 


1.0000(1) 


1.31(1) 


3.75 


4.0000(1) 


1.24(1) 


3.0000(1) 


1.24(1) 


2.0000(1) 


1.24(1) 


1.0000(1) 


1.25(2) 


4.00 


4.0000(1) 


1.20(2) 


3.0000(1) 


1.19(1) 


2.0000(1) 


1.19(1) 


1.0000(1) 


1.20(2) 


4.25 


4.0000(1) 


1.13(2) 


3.0000(1) 


1.13(1) 


2.0000(1) 


1.13(1) 


1.0000(1) 


1.14(1) 


4.50 


4.0000(1) 


1.10(2) 


3.0000(1) 


1.10(1) 


2.0000(1) 


1.10(2) 


1.0000(1) 


1.11(2) 



11(b) Critical values for energy Ec and exponent ve for different a and disorder strength W 



a 


W=3.0 


W=2.0 


W=1.0 


E c 


ve 


E c 


ve 


E c 


ve 


2.25 


0.5007(4) 


2.15(4) 


1.0004(3) 


2.14(3) 


1.5003(4) 


2.12(4) 


2.50 


0.5003(2) 


1.84(2) 


1.0002(2) 


1.85(2) 


1.5001(2) 


1.84(3) 


2.75 


0.5001(1) 


1.66(1) 


1.0001(1) 


1.65(1) 


1.5000(1) 


1.66(1) 


3.00 


0.5001(1) 


1.51(1) 


1.0001(1) 


1.51(1) 


1.5000(1) 


1.52(1) 


3.25 


0.5000(1) 


1.40(1) 


1.0000(1) 


1.41(1) 


1.4999(1) 


1.43(2) 


3.50 


0.5000(1) 


1.30(1) 


1.0000(1) 


1.31(1) 


1.5000(1) 


1.32(1) 


3.75 


0.5000(1) 


1.24(1) 


1.0000(1) 


1.24(1) 


1.5000(1) 


1.25(1) 


4.00 


0.5000(1) 


1.19(1) 


1.0000(1) 


1.19(1) 


1.5000(1) 


1.19(2) 


4.25 


0.5000(1) 


1.14(1) 


1.0000(1) 


1.14(1) 


1.5000(1) 


1.15(1) 


4.50 


0.5000(1) 


1.10(1) 


1.0000(1) 


1.10(2) 


1.5000(1) 


1.11(2) 



TABLE II: The values of critical exponent v and critical points Wc and Ec got from the fitting A with Eq. (10) for a > 2. The 
errors are obtained from bootstrap method with 95% confidence. 
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FIG. 6: The a-dependence of the averaged critical exponent 
v (for different values of W c and E c ) and the fitting curve 
from Eq. (11) with Ai = 1.4 and Ai = 1. Inset: ln(^ — 1) as 
a function of a — 2, where the straight slope is set to —1. 



ical parameters. There is a finite range of extended 
eigenstates, and the effective bandwidth decreases lin- 
early with increasing W for a > 2. The positions of the 
mobility edges separating localized and extended states 
depend only on the disorder strength whenever a > 2. 
Using the finite-size scaling analysis, we find the critical 
exponent v = 1 + 1.4exp(2 — a). In particular v = 2.4 
when a = 2 and v = 1 when a — > oo. Compared to the 
3D uncorrected Anderson model, the positions of the 
mobility edges are apparently different. And the critical 
exponent v is no longer equal to 1.58 17,31 . So, the transi- 
tion in the ID disordered system with correlation is in a 
new universal class. The existence of mobility edge and 
a dependence of the critical exponent indicate that the 
nature of a correlated disordered system is somehow be- 
tween that of the pure random (without any correlation) 
and pseudo-random (quasi-periodic) systems. Our work 
sheds some light on the metal-insulator transitions. 



IV. CONCLUSION 
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